Giant vortex and Skyrmion in a rotating two-species Bose-Einstein condensate 
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Numerical simulations are performed for a rotating two-species Bose condensate confined by a 
harmonic potential. The particle numbers of each species are unequal. When the rotational speed 
exceeds a critical value, the majority species reside in the center of the potential while the minority 
species is pushed out to the outskirts, forming a giant vortex hole to contain the majority species. 
A novel annular Skyrmion forms at the interface of the two species. 
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INTRODUCTION 

Quantum coherence enables intriguing phenomena in 
Bose-Einstein condensates (BEC) as quantized vorticity, 
which has been verified in many experiments. When a 
trapped BEC is driven to rotate, singly quantized vor- 
tices form. Faster rotation generates more vortices which 
finally condense into a lattice. Other methods, such as 
the phase imprint technique, are also used to create vor- 
tex in BECs. D.R. Scherer et al[l] carried out an exper- 
iment to implement vortices by the interference of three 
independent trapped BECs. 

There are many efforts to create a multiply quantized 
vortex or giant vortex in the BECs. However, a giant 
vortex is energetically disfavored in a harmonic trap and 
is not expected to persist if created in a rotating con- 
densate. Some authors tried to overcome this dissocia- 
tion instability by employing an external repulsive pin- 
ning potential [2[ or a steeper trap that allows rotation 
at frequency exceeding the centrifugal limit of the har- 
monic trap. [3, 0, [H, @, uj Multiply quantized vortices are 
also created topologically by deploying the spin degrees 
of freedom in optical traps. [8[ In a recent experiment, P. 
Engels et al shone a near resonant laser beam through a 
rapidly rotating harmonically trapped BEC to produce a 
density hole encircled by a high number of vorticity. [9|, LL0[ 

In this work, we produce a giant vortex in a harmonic 
trap by introducing a second species of BEC. There is a 
population imbalance (iVi > A^) in the two- species BEC 
and the inter-species scattering length is larger than the 
intra-species length. We consider a quasi-2D repulsive 
system by numerically integrating the Gross-Pitaevskii 
(GP) equation. S. Bargi et al[l6[ have studied the same 
setup by numerical diagonlization and variational mean- 
field approxiamtion. In the limit of Ni ^> A^, the density 
profile of the majority species is approximately the same 
as if the minority species is absent. On the other hand, 
the minority species feels a trap formed by the external 
magnetic trap plus the hump in the center which is gener- 



ated by the majority species. Hence the effective poten- 
tial acted on the second species behaviors like a 'Mexican 
hat'. The minority species is pushed to the rims, which 
makes for the formation of a hole. When this system 
is set to rotate, a series of singly quantized vortices are 
created. Large centrifugal force induced by faster rota- 
tion subsequently generates a giant hole in the minor- 
ity species to contain the majority species. A number 
of phase defects congregated in the central area reveals 
that this giant hole is a multiply quantized vorticity. This 
becomes possible because of the low density in the hole 
costing little energy for vortices to merge into the giant 
vortex. By describing the system in terms of pseudo- 
spin density parameters, it is found that associated to 
the giant vortex is an annular Skyrmion. Although there 
are many works to deal with Skyrmions resulted from 
singly quantized vortices, [III [H, |H, EH the present 
multiply quantized vortex creates a Skyrmion with novel 
and more complex topological structure. The topological 
charge density exhibits a ring distribution. 



FORMULISM 

In weak interaction limit, quantum dynamics of the 
BEC is governed by the nonlinear Gross-Pitaevskii equa- 
tion. We consider a 2D system that subjects to rotation 
Q = Qz. In a frame rotating with frequency ft around 
the z axis, the dynamics of a condensate wavefunction 
if = (^i,^2) T is described by the following normalized 
equation 



V(f)+J2 9ijHj\ 2 -^L z }A, (1) 
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where gij,(i,j — 1,2) are the nonlinear coupling con- 
stants which are expressed by — AnNjaf^ with afj 
the inter (i ^ j) or intra (i = j) species s-wave scat- 
tering lengthes. V(f) = \uo 2 (x 2 + y 2 ) is the harmonic 
trapping well. We assume that the two species of BEC 
are trapped by the same external potential with unequal 
particle numbers N\ > N2 and all the nonlinear coupling 
constants are positive. 
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We first analyze the problem in the limit of N\ ^> 
N 2 . The distribution of the majority species is hardly 
affected by minority species. Hence we can approxi- 
mate its profile with the Thomas-Fermi (TF) distribu- 
tion, m(f) = |^i(r)| 2 = [/ii - \(uo 2 - n 2 )r 2 ]/ gil for 
n\(r) > and n\(r) = otherwise. is the chem- 
ical potential. The term associated with the rotation 
frequency ft is the centrifugal potential. We consider 
the case for Q < uj. The radius R\ of the conden- 
sate is evaluated by finding where the density goes to 
zero. If iii > then there is a single solution with 
Ri = Y // 2/ii/(o; 2 — ^ 2 ). The chemical potential is deter- 
mined by normalization condition J ni(r)df = 1, which 
gives rise to \i\ = y gn(oJ 2 — Q 2 )/7T. Because of the 
strongly repulsive interspecies coupling, the minority is 
pushed out to the rim which feels an effective potential 
produced by the external trap combining with the repul- 
sive hump of the majority, V e ff(r) = V(r) + #21^1(^)5 
which looks like a 'Mexican hat'. The minority conden- 
sate prefers to reside in annular notches of the effective 
potential and subsequently forms a huge hole in the cen- 
tral area to contain the majority condensate. 

The Thomas-Fermi profile of the minority condensate 



n 2 (f) = \Ur)V = (M2 - V eff (r) + -fiV)/ ff22) (2) 

for ri2 (r) > and 712(f) = otherwise. The outer and 
inner radii of the ring-shape condensate are 
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The corresponding chemical potential is given by 
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which requires #21 > #11 • In this case, the two species of 
BEC is imiscible. 

Suppose there is a vortex in the center with circula- 
tion T = 27rm. The quantization condition requires the 
vorticity of integer values m = 0, 1, 2, • • •. The energy 
of the minority condensate in the TF approximation is 
expressed as, 



E 



J df\Veff 



TTL 1 

+ ^2 + o#22n 2 (r)]n 2 (r). (5) 



In the above analysis we have omitted the contribution 
from singly quantized vortices which may form a vortex 
lattice. The energy in a frame rotating with angular ve- 
locity Q is related to that in the laboratory frame by 
E'(m) = E(m) — mVt. For a given rotation frequency 
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FIG. 1: The rotational frequency dependence of total vorticity 
from formula [6]. 



£7, this energy is to be minimized with respect to the 
parameter m to obtain the vorticity of the giant vortex. 
By treating the parameter m as continuous variable, the 
final result is 



m 



922^ 



MM2ln(f)-^Miln(f) + i(^ 
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(6) 



Figure 1 depicts the vorticity dependence on the rota- 
tional speed Q. It should be noted that the total vorticity 
also includes the singly quantized vortices. Hence Eq.[6] 
is a schematic formula for the relation between the vor- 
ticity and the rotational speed. 



NUMERICAL SIMULATIONS 

The numerical simulations is carried out by solving the 
norm-preserving imaginary time propagation of Eq. [I] . 
The initial wavefunction adopts the Thomas-Fermi ap- 
proximation. We adopt the time-splitting Fourier pseu- 
dospectral method developed by Bao et al[17] to com- 
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FIG. 2: Density profiles of the two- species BEC for rota- 
tional speed Q/u; = 0.85. The ratio of particle numbers is 
N1/N2 = 4. In the minority species there forms a ring of 
singly quantized vortices plus a huge hole in the center. 

pute the partial differential equation [I] in a region 
0, y) e [-8, 8] x [8, 8] with a refined grid of 80 x 80 
nodes, which is sufficient to achieve grid independence. 
The particle number ratio is fixed at N1/N2 = 2 and the 
nonlinear coupling parameters are chosen as gu = 100, 
#22 = 550, 012 = 200, and g 21 = 400. 

Figure 2 shows the density profile of both species of 
condensates. As we have expected, the majority con- 
densate resides in the central area of the harmonic well 
while pushes the minority condensate to the outskirts, 
enabling the latter to acquire more angular momentum 
to circulate around it. In the mean time, the majority 
condensate becomes more compact by the inward force 
implemented by the minority condensate. As the rotat- 
ing speed increases, a lot of phase defects of the minority 
condensate enter the center area of the harmonic trap 
and merge into a giant vortex. Here the giant vortex 
implies that several phase defects are contained in a sin- 
gle density hole. Figure 3 exhibits that the circulating 
movement of the minority condensate forms a giant vor- 
tex with winding number m = 8. In addition, a sequence 
of singly quantized vortex forms along the ring-shape mi- 
nority condensate. 




FIG. 3: Phase plots for (a) the majority species and (b) the 
minority species. The parameters are the same as in Fig.l. 
Panel (b) indicates that the huge hole is a multiply quantized 
vortex with a winding number m = 8. 

To further reveal the explicit structure of the giant 
vortex, we introduce a normalized complex valued spinor 
X — [Xi(r)j X2(f)Y by decomposing the wavefunction 
ipi = yPr{f)Xi{^)^ where pr is the total density. They 
satisfy |xi| 2 + \X2\ 2 = 1- The pseudo-spin density is de- 
fined as S = X^(r)<?x(r)> wnere & is the Pauli matrix. 



Obviously, the modulus of the total spin is | aS' | = 1. The 
pseudo-spin texture corresponding to the giant vortex is 
shown in Fig. 4. It is shown that the S x or S y exhibits an 
m-fold symmetrical annular necklace of radius i?, which 
is clearly related to the quantized circulation of the giant 
vortex. As S x = 2|xi I IX2 1 cos(#i — #2), where 6i(f) are 
the phases of the condensate wavefunctions, and from 
Fig. 3 the phase of the majority condensate is approxi- 
mately constant, we know that S x has m-fold symmetry. 
However, if the rotation frequency Q increases further, 
vortices are created in the majority condensate and 0\ 
can not be viewed as constant, then the m-fold symmetry 
should be broken. Fig. 4(b) is the S z distribution which 
indicates that the pseudo-spin points up at the center 
and points down on the edge which results from \i ~ 
for r < R and xi ~ for r > R. At the interface of the 
two condensates, the system twists its pseudo-spinor or- 
der parameter, as shown in the projected vectorial plot of 
(S x ,S y ) in Fig. 4(c). It reveals a rather complex pseudo- 
spin texture which may be called a giant Skyrmion. 




FIG. 4: Pseudo-spin density distribution for (a) S x , (b) S z . 
Panel (c) is the vectorial plot of (S x ,S y ), which exhibits an 
annular Skyrmion structure, (d) is an amplified local part of 
(c). 

The pseudo-spin texture in Fig. 4 reveals a concrete 
/c-fold symmetry. At first sight, it seems curious that 
k = m— 1 instead of k = m. It is because the pseudo-spin 
contributes is an additional rotational angle A(/> = 2ir/k 
that relates to the spatial rotation. Taking the spins 
to rotating about the z axis, the symmetry is formally 
described by 

5 i (i?(27r/fe)f) = 7^(27T + A0)S,-(f), (7) 
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where 



R{4>) 



COS <p 
— sin ( 



sin c 

cos ( 



(8) 



represents a spatial rotation about the z axis by an an- 
gle 0, and lZij{(j)) = x, y) represents a pseudo-spin 
rotation about the z axis by an angle (j). The projected 
pseudo-spin vector increases an angle of 2i\m when it 
runs along a route that encloses the origin. 

The left panel of Fig. 5 shows the effective velocity 
field that is defined as 



Veff(r) = (3i(r)+n(r))/pT(r), 



(9) 



with ji 57 00* V - fa V ip*), 1, 2) the partial 
current density. Besides a sequence of singly quantized 
vortices, the whole minority condensate circulates around 
the majority condensate with a multiply quantized circu- 
lation r = 8 x 2tt. Contrary to the conventional vortex in 
a single species condensate, \v e ff\ vanishes at the center 
which makes a coreless vortex without a density dip in 
the total density. The largest curl of the velocity field 
lies on a ring of finite radius. 
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To simplify the problem, we omit the chain of singly 
quantized vortices and focus our attention on the giant 
vortex. The configuration satisfies the boundary condi- 
tion /3(0) = and /3(oo) = 7r, which is referred as the 
Anderson- Toulouse vortex. The phase of majority con- 
densate is constant which can be set as 0\ = 0. The 
phase of minority condensate is independent of radius r 
and can be approximately written as 62 = m<\>. Hence 
the spatial rotation of an angle 2i\jk corresponding to a 
relative phase difference A6 = 2i\mjk. 

The pseudo-spin vector S can be expressed in the polar 
coordinates 



5 = Wl-[S z (r)] 2 cos(m<p), y/l - [S z (r)] 2 sin(m^), S z }. 

(12) 

The topological charge density of Eq. [10] has the form 



q{f) 



m dS z (r) 



(13) 



4nr dr 

which leads to the total charge of the giant Skyrmion as 



m 



d*rq{r) = -[5 z (oo) - S z (0)] = m. (14) 



In summary, we have created a giant vortex in a two- 
species Bose condensate with unequal particle numbers. 
The corresponding pseudo-spin texture shows a m-fold 
symmetry and the topological charge is m. 
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FIG. 5: (Color online) (a) The vectorial plot of the effective 
velocity v e ff defined by Eq.[9]; (b) The topological charge 
density q(r) distribution. 



The topological charge density is defined as 
1 



q{f) 



8tt 



-e lJ S-diS x dnS. 



(io) [3; 



The topological charge density q(f) characterizes the spa- 
tial distribution of the Skyrmion. The total topologi- 
cal charge or the Pontryagian index Q = J drq(r) is an 



invariant. [18l UJj From the right panel of Fig. 5, one 
finds that the singly quantized vortices have the Dirac 5- 
function topological charge density. On the other hand, 
the giant Skyrmion has its charge uniformly distribute 
on a ring where the two imiscable condensates overlap. 



DISCUSSIONS 

The additional rotation angle A<p can be accounted for 
by parameterizing the wavefunction as 



Xi 

X2 



cos/3(r)/2 \ 
! sin/3(r)/2 J ' 



(11) 
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